Thermodynamics and Fluctuation Theorems for a Strongly Coupled Open Quantum 

System: An Exactly Solvable Case 
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We illustrate recent results concerning the validity of the work fluctuation theorem in open quan- 
tum systems [M. Campisi, P. Talkner, and P. Hanggi, Phys. Rev. Lett. 102, 210401 (2009)], 
by applying them to a solvable model of an open quantum system. The central role played by 
the thermodynamic partition function of the open quantum system, - a two level fluctuator with a 
strong quantum nondemolition coupling to a harmonic oscillator -, is elucidated. The corresponding 
quantum Hamiltonian of mean force is evaluated explicitly. We study the thermodynamic entropy 
and the corresponding specific heat of this open system as a function of temperature and coupling 
strength and show that both may assume negative values at nonzero low temperatures. 

PACS numbers: 05.30.-d, 05.70.Ln, 05.70.-a, 65.40.gd, 65.40.Ba, 



I. INTRODUCTION 

Exact results about nonequilibrium fluctuations in 
nanosystems, such as the Jarzynski equality [l| and the 
Tasaki-Crooks fluctuation theorem have recently at- 
tracted a great deal of attention in the burgeoning field of 
nonlinear fluctuation relations. These results were first 
derived for classical systems and later for quantum sys- 
tems that are either isolated or weakly coupled to their 
environment @, 0, i, @, 0, 0, i, 0, E(- However, often 
the interaction with the environment does play an impor- 
tant role which cannot be neglected in real experimental 
situations. For this reason more attention has been re- 
cently devoted to the strong coupling regime, both classi- 
cally [l2| , and quantum mechanically [l3j . In this regime, 
the driven system of interest (with Hamiltonian Hs(t)), 
strongly couples to a bath (H B ), via a non negligible 
interaction term Hsb- 



H(t)=H s {t)+H SB +H B 



U) 



For the applicability of work and fluctuation theorems, 
the correct choice of the statistical mechanical descrip- 
tion of an open quantum system in terms of the proper 
thermodynamic partition function, is of decisive impor- 
tance. The bare system partition sum 



Qs(t) = Ti s e 



-0Hs(t) 



(2) 



clearly fails to account for the effects of the environment 
on the system. It rather is the open system thermody- 
namic partition function 



Zsit) 



Y(t) 



(3) 



denotes the total system partition function; i.e. 

Y(t) = ^re-PiHsi^+HsB+HB) 
and Z B the bare bath partition function 
Z B = Tr B e-^ B 



(4) 



(5) 



The time t merely specifies the values of the external 
parameters as they occur in the course of the driving 
protocol at the time t. The symbols Trg, Tr B , Tr denote 
traces over system, bath, and total system respectively. 
The symbol (3 = (k B T)~ l indicates the inverse thermal 
energy, with k B Boltzmann constant and T the tempera- 
ture. This temperature is provided via vanishingly small 
weak contact with a large (super)-bath, which allows for 
a statistical mechanical treatment. 

The adoption of the thermodynamic partition func- 
tion Zsit), and the corresponding free energy Fs(t) = 
— /3 _1 \aZs{t) allows to obtain the Jarzynski equality 



-0AF s 



(6) 



valid irrespectively of the coupling strength [l3| . 

In the following we exemplify this result by applying 
it to a simple model Hamiltonian of an open quantum 
system, Sec. HH Mil We next illustrate the equilibrium 
thermodynamics of that open system, by computing its 
Hamiltonian of mean force, its entropy and specific heat, 
see Sec. IIVI Remarks and conclusions are drawn in Sec. 

El 



II. A TWO LEVEL 
FLUCTUATOR-OSCILLATOR MODEL 

We consider the following Hamiltonian describing a 
two level system and a harmonic oscillator interacting 
with each other: 



H(t) = ^<t 2 ( 



>l B +l s ®tt 
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(7) 
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Here u z is a Pauli matrix of the two level system, a) 
and a are raising and lowering operators of the harmonic 
oscillator, e(t), tt, \ are the two level system energy spac- 
ing, the oscillator energy quantum and the coupling en- 
ergy, respectively. The parameter x can assume positive 
and negative values whereas e(t) and f2 are strictly pos- 
itive. The oscillator energy quantum f2 is related to the 
oscillator frequency uj, via Planck's constant £1 = Huj. 
We consider the two level system (also referred to as 
the qubit throughout the text) as our system of inter- 
est (Hs(t) = e(t)a z /2 ® Is), and the oscillator as our 
stylized, "minimal" bath; i.e. Hb = fl(a^a + 1/2) $5 lg. 
The operators Is and 1b denote the identity operators 
acting on the system and bath Hilbert spaces, respec- 
tively. We require |%| < ft, which ensures that the total 
Hamiltonian is bounded from belowguaranteeing stabil- 
ity of the total system (from Eq. ([9|), £1 + \s must be 
positive in order that the smallest eigenvalue be finite). 
The two level system energy spacing e is assumed to de- 
pend on time according to some pre-specified protocol. 
This model Hamiltonian has the peculiarity that the in- 
teraction Hamiltonian commutes with both system and 
bath Hamiltonians, implying a so called quantum nondc- 
molition coupling: 

[H s (t) ® i B , H SB ] = [Is ® H B ,H SB } = . (8) 

The time-instantaneous energy eigenvalues assume the 
form 



£(*), 

2 











(raH 


4) 


+ X s ( n - 


-I) 



(9) 



s = ±1, n = 0,1,2,... . We remark that the corre- 
sponding instantaneous eigenstates |n, s) do not depend 
on time. 

The partition function Y(t) of the total system be- 
comes, with Eq. (fj| 

V/^-M = ?+ (t) +? _(t) (10) 



where 



9±W 



e -/3n/2 eT/ 3(e(t)+ x )/2 

1 _ e -P(n± x ) 



(11) 



The bare bath partition function Zb is, with Eq. ([5]): 



,-/3fi(n+l/2) _ 



2sinh(/?rj/2) 



(12) 



Then the thermodynamic partition function Zs(t) of the 
open system becomes, according to Eq. ([3]): 

Z s (t) = 2(q+(t) +q-(t))smh(f3n/2) . (13) 

Note that the open system thermodynamic partition 
function differs substantially from the bare system parti- 
tion sum Qs{t), which reads, with Eq. ([2]): 

Qs(t) = J2 e- pe(t)s/2 = 2 cosh(/?e(t)/2) . (14) 



In particular the thermodynamic partition function con- 
sistently accounts for the presence of the oscillator and 
the interaction, as it depends on £1 and x-i whereas the 
partition sum Qs(t) does not. 



III. WORK AND FLUCTUATION THEOREMS 

For a prescribed protocol of the two level spacing e(t) 
to < t < tf, the work w performed on the two level 
system is distributed according to the probability density 



ingt 
S 



function Pt fl t (w), given by 

oo 

Pt f ,t (w)= H S ( w ~ ( E ™Atf) -E n , s (ta))) 

e -0E n , s (t o ) 
Y(t ) 

(15) 



m,n— r,s— ±1 



x P(m, r\n, s)- 



wherc 8(x) denotes the Dirac delta function and 
P(m,r\n, s) is the transition probability to jump from 
the eigenstate \n, s) of the total Hamiltonian at time to 
to the eigenstate \m, r) at time tf: 



P(m,r\n,s) = \(m,r\U tf ,to\n,s)\' 



(16) 



with Ut f ,t = Texp(— i J dtH(t)/h) denoting the time 
evolution operator. The model Hamiltonian in Eq. I|7|). 
commutes with itself at different times, so the time or- 
dered exponential reduces to an ordinary exponential 



exp 



, £ (*) - 

dt—a z 

2 z 



at a + - | (tf-to) 



(17) 



By inserting this expression into Eq. (|16[) . one sees that 
no transition takes place 



P(m,r\n,s) = S mn 6 r 



(18) 



with 5 m ,n denoting the Kronecker symbol. This is of 
course to be expected since an interaction that commutes 
with the free evolution does not cause any transition. 
Thus, for the work probability density one obtains: 

Pt f ,t ( w ) = 

q+(t )S(w - Ae/2) + g-(t )6(w + Ae/2) (19) 
q+{t ) +q-(t ) 

where Ae = e(i/) — e(to). By exchanging to with tf one 
obtains the backward pdf of work pt ,tf (w), correspond- 
ing to the backward protocol e(t) = e(tf + to — t). After 



3 



some calculations one obtains the following expression for 
their ratios: 



o 0w 



cosh(/3 



e(«/)+X 



-0CI 



cosh(/3 



-(tf)-x 



cosh(/3^±^) 



e~f m cosh(/?2(£^) 



Ptj,t (w) 
Vt Q ,t f {-w) 

"(20) 

where we recognize that the ratio on the right hand side is 
equal to Zs(tf)/Zs{to), as predicted by the work fluctu- 
ation theorem for arbitrary open quantum systems [HI]. 
Using (JTHJ) we also obtain the following expression for the 
Jarzynski exponentiated work: 



-0W\ 



cosh(/3^±^) 



cosh(/3 



g(«/)-x 



cosh(/3 



e(*o)+X 



cosh(/3 



l(*q)-2 



(21) 
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n_ 

~2e 



that is, 



(e-*°) 



^s(*o) 



(22) 



as predicted by the Jarzynski equality for arbitrary open 
quantum systems in Eq. ([6]) [1 31 ] - 

By comparison of Eqs. (fT4")l and (|2"Tj) one observes that: 



Qs(t f ) 

Qs(t y 



(23) 



This result is in contrast to recent claims reported by 
Teifel and Mahler [28|, according to which the averaged 
exponentiated work should be identical to the ratio of 
partition sums Qs independently of coupling strength, 
provided the interaction commutes with both system and 
bath Hamiltonians as it is the case with the present study. 



IV. EQUILIBRIUM THERMODYNAMICS 

We turn now to the study of the equilibrium thermody- 
namics of the open two level system. This means that we 
now keep e fixed and study the time independent Hamil- 
tonian 
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FIG. 1: Dimensionless difference between renormalized 
qubit's energy spacing and original qubit's energy spacing, 
A/s (solid line, Eqs. 1)28 130 and global dimensionless shift 
of the energy spectrum, j/e (dashed line, Eq. 1)29 J l). as func- 
tions of (a) the dimensionless coupling strength \/e (top 
panel), (b) dimensionless oscillator's energy quantum fi/e 
(bottom panel). The graphs correspond to a temperature 
of T = 50mK, and the experimental values employed in Ref. 
[33 e/2-Kh = 6.9 GHz, (a) Q,/2nh = 5.7GHz (top panel), (b) 
x/nh = -17MHz (bottom panel). 



(24) 



The Hamiltonian of Mean Force 



A fundamental quantity that is closely related to the 
open system partition function Zg is the quantum Hamil- 
tonian of mean force [l3[ 



1 Trvr-P^s+HsB+HB) 

H* := -A m ™ . 

a z B 



(25) 



It generalizes the potential of mean force commonly em- 
ployed in reaction rate theory (29| and implicit solvent 



models [30|. The Hamiltonian of mean force is the ef- 
fective Hamiltonian that describes the open system at 
equilibrium with the environment according to the equa- 
tion: 



Z^e 



-I3H 



(26) 



It hence determines the reduced density matrix of the 
open system, ps, in thermal equilibrium according to 
ps = Zg l e~P H . The calculation of H* in general is 
a difficult task. However for the model Hamiltonian in 
Eq. (|24|) . the calculation is straightforward and leads to 



h* 



(27) 



4 



where 



£1/8 



e + x + 77 artanh 

P 



1 



~P n sinh(/?x) 
e -pn CO sh(/3x) 



is the renormalized level spacing, and 



7 



Am 

2/3 



1 - 2e- 0n cosh(/3 x ) - 
(1-e-/ 30 ) 2 



5 -2/3n 



(28) 



(29) 



specihes a global shift of the spectrum. In obtaining Eq. 
(|27)l we used the identity e aa * = cosh(a)l,s + sinh(a)i7 2 . 
When x ~ * 0, the renormalized spacing tends to the orig- 
inal spacing e, and the offset 7 vanishes, so that H* tends 
to the bare system Hamiltonian H$, as expected. Fig. 
QJa) displays 7 as well as the amount of renormalization 

A := s* - s (30) 

which is independent of the bare spacing e, as functions 
of the coupling strength x, for |x| < f2. These quantities 
are displayed in non-dimensional units where energies are 
rescaled by e. As \x\/ £ approaches the stability limit 
■iSi/e, A/e and 7/e diverge, while they vanish as the 
coupling x/ £ approaches zero. From Eq. (|28[) . we note 
that, given certain values of the spacing e and of Q, there 
exists a value of x f° r which the renormalized energy 
spacing e* vanishes, meaning that an effective degeneracy 
of the qubit is induced by the presence of the oscillator. 
In Fig. DJb), A/e and 7/e are plotted as functions of 
Q/e, for fixed (3 and Xi an d for fi/e > |x|/e- 

The graphs in Fig. [1] correspond to values of f2 and x 
that match the regime of values used in an experimental 
implementation of the model Hamiltonian in Eq. (|24p 
with superconducting circuits, as it has been recently re- 
ported [311 ]. In that experiment |x| is about two orders of 
magnitudes smaller than Q, thus the leading corrections 
to the energy spacing are of first order and those of the 
shift 7 are at most of second order in x/ £ - 

At low temperatures we find the following limiting re- 
sults: 



lim 7 
lim £* 





e + x 



gk 

lim — — r£ 

/3^oo d(i k 



—r;£* = k = 1,2,3.. 



(31) 
(32) 

(33) 



From the previous two equations we deduce that the de- 
generacy of the spectrum occurs at T = for the special 
value x = ~ e - I 11 the following we will come back to 
the effect of this degeneracy on the system's entropy and 
specific heat. 



B. Thermodynamic Entropy 

From the partition function, Zg, one obtains the free 
energy: 

F s = -k B T\nZ s = -{l/P)\nZ s , (34) 



■Sl/e 




FIG. 2: Contour plot of dimensionless entropy Ss/kB, Eq. 
(|35[) . as a function of dimensionless temperature fcsT/e = 
l/(/3e), and dimensionless interaction strength x/ £ y f° r ( a ) 
Q/e = 3 (top panel), (b) Q/e = 1/3 (bottom panel). The 
entropy is nowehere negative in panel (a). In panel (b), it 
assumes negative values in the region, labelled as Ss < 0, 
enclosed by the level line Ss — (thick light blue line). 



and the entropy: 
Ss 



dFs R2 dF s 
= ksp 



dT 



Oil 



(35) 



In Fig. HJa) the entropy following from Eq. (JT3J) with 
e(t) = e, is displayed as a function of dimensionless tem- 
perature ksT/e = \/{f3e) and dimensionless coupling 



5 



strength, \l £ for a fixed value of rescaled oscillator en- 
ergy quantum fi/e, larger than 1. As \/e approaches 
the instability values ±f2/e the entropy diverges. For all 
values of x/ £ i the entropy vanishes at zero temperature 
in agreement with the third law. An exception is at the 
special case \l £ = — 1) where the ground state assumes 
a finite degeneracy. Put differently, for x/ e = ~ lj the 
zero temperature entropy is no longer zero but assumes 
the finite positive value fc_B In 2 [32|, [H, H3, [35[ . This 
Ub In 2 term is a consequence of the fact that in the limit 
of zero temperature and for x = ~ £ i the effective spac- 
ing £* and all its higher order derivatives with respect 
to temperature vanish (see Eq. (|32l33p ). At finite tem- 
peratures there are values of x f° r which the spacing e* 
vanishes, however these do not coincide with the values 
of x f° r which the entropy is fcs In 2, since then the first 
derivative of e* with respect to (3 does not vanish and 
consequently yields a contribution to the entropy. 

Fig. [D(b) depicts the entropy as a function of dimen- 
sionless temperature ksT/e = l/(/?£) and dimensionless 
coupling strength, xl £ i f° r a fixed value of rescaled os- 
cillator energy quantum Sl/e, smaller than 1. The most 
prominent feature in this case Q/e < 1 is the appear- 
ance of a region of negative entropy for small values of 
ksT and negative coupling (the region labeled as Ss < 
in Fi g. Ef b)). Interestingly, the experiment reported in 
Ref. [31| lies in this region SI < e, x < where the en- 
tropy may become negative. For the parameter values 
reported therein, a negative entropy is expected below 
~ 22 mK. From Fig. Hlb) we notice that, for positive 
X, the entropy vanishes at absolute zero temperature, in 
accordance to the third law of thermodynamics [32j , and 
reaches a plateau at high temperatures, without becom- 
ing negative. For negative x, it vanishes as well at zero 
temperature, however with increasing low temperatures, 
entropy first decreases until it reaches a negative mini- 
mum, and then increases until it approaches a positive 
plateau-value at high temperatures. This behaviour is 
further illustrated in Fig. [31 

Independently of the sign of \i a t high temperature 
the entropy reaches the asymptotic value: 



lim Ss 

13^0 



kn In 



2Q 2 



W-x' 



(36) 



which notably does not depend on the spacing e. For 
X = 0, the high temperature entropy in Eq. (f36|) be- 
comes equal to ks In 2, reflecting the fact that spin up 
and spin down states become equally populated at infi- 
nite temperature. For x 7^ 0, it assumes values larger 
than A:sln2, and diverges for |x| approaching £1. 



C. Specific heat 



From the entropy (|35|) one obtains the the specific heat 
of the open two level system: 



C s 



T 



dS s 
dT 



-0 



dS s 

d(3 



1 

0.75 



-£ 0.5 



0.25 














' — X/e = 


-1/6 - 




---- X/e = 


1/6 . 











0.5 



1.5 



k B T/e 



(37) 



FIG. 3: Dimensionless entropy Ss/ks as a function of dimen- 
sionless temperature ksT/e — l/(/3e), for Q/e = 1/3, and 
X/e = —1/6 (solid line), \/e — 1/6 (dashed line). The dotted 
line is the asymptotic value calculated from Eq. (1361 1, 



Figs. [Ha) and3Jb) represent the specific heat as a func- 
tion of dimensionless temperature ksT/e = l/(f3e) and 
dimensionless coupling strength, x/ £ f° r fixed values of 
rescaled oscillator energy quantum il/e. 

In Fig. HJa) fl is larger than e. For values \x/ £ \ < 
fi/e, the specific heat vanishes at zero temperature. With 
growing temperatures, it first increases, then reaches a 
maximum and finally goes again to zero. The maximum 
occurs at decreasing temperatures as x/ e ~ ¥ ~ 1- As x/ £ 
approaches —I we also see the appearance of a second 
maximum for larger values of fc^T. These features are 
further illustrated in Fig. EJa). 

The specific heat landscape changes drastically for 
ft < e, Fig. 0Jb). The most relevant feature in this 
parameter range is the appearance of a region of nega- 
tive specific heat at low temperatures and negative xl £ 
(the region labelled as Cs < in Fig. Sfb)). From 
Fig. Bib) we observe that, for < xl £ < ^V e the spe- 
cific heat starts from zero at zero temperature, reaches a 
maximum and goes to zero again at high temperatures. 
For — Vt/e < x/e < the specific heat starts at T = 
from zero, reaches a negative minimum with increasing 
temperature, then a positive maximum and finally goes 
to zero at high temperature. These features are further 
illustrated in Fig. [^b). From Fig. \5j(b) we also no- 
tice that the curves corresponding to positive x all cross 
each other within a very small temperature range around 
ksT/e ~ 0.21 for the given value f2/e = 1/3. Indeed, for 
ksT/e = 0.21, fi/e = 1/3, and xl £ ranging from to 
f2/e, the specific heat is almost constant (with variations 
within 5 % of its value). An analogous situation hap- 
pens also for other values of fl/e < 1, showing that in 
this regime one should expect the presence of a narrow 
temperature range for which the specific heat is not very 
sensitive to changes in coupling strength Xi a s long as 
this remains positive. Regardless of whether is larger 
or smaller than e, the specific heat approaches a unique 
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-n/e 



Q./E 




-n/e 



FIG. 4: Contour plot of dimensionless specific heat Cs/k,B, 
Eq. (|37[) . as a function of dimensionless temperature ksT/e — 
l/(/3e), and dimensionless interaction strength x/ e > f° r ( a ) 
fi/e = 3 (top panel), (b) fi/e = 1/3 (bottom panel). The 
specific heat is nowehere negative in panel (a). In panel (b), 
it assumes negative values in the region, labelled as Cs < 0, 
enclosed by the level line Cs = (thick light blue line). 



functional form in the limits as x approaches ±fi. This 
limiting function can be calculated analytically: 



lim Cs = ki 



2k B T 
ft 



sinh 



2k B T 



(38) 



The fact that it does not tend to zero at zero temperature 
is not in contrast with the third law of thermodynamics, 
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0.5 0.75 
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FIG. 5: Dimensionless specific heat Cs/ks as a function of 
rescaled temperature fcsT/e = l/(/3e), for various values of 
X/e and (a) fi/e = 3 (top panel), (b) fi/s = 1/3 (bottom 
panel). 



since for x = the system is no longer stable. 

Fig. [S] depicts the behavior of the specific heat as 
x/e — ► ±Q/e for fi/e < 1. The convergence to the 
limiting function in Eq. |38|) as x/e — > +fl/e is quite 
fast compared to the much slower convergence in the 
other limit x/ £ ~ * — f2/ff- The approach to the limit 
is qualitatively very different in the two cases. In both 
cases the specific heat vanishes at low temperature and 
it approaches the limiting curve in Eq. (f3"B")) for large 
temperatures. However in the case x/ £ ~ ~^V £ ; the 
specific heat displays a drastic peak at intermediate tem- 
peratures. As x/e approaches —Q/e, the peak becomes 
increasingly pronounced while getting closer to the ori- 
gin of the temperature axis. In the limit x/ e ~ * ~ ^/e, 
eventually the peak becomes a delta singularity at zero 
temperature. This singularity contributes with a finite 
term A to the total heat Q = J CsdT, which, in analogy 
with first order phase transitions, can be interpreted as 
a latent heat. 
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\ ---- x/e = -n/e + W- 7 

\ x/e = -n/e + 10- 10 

\ x/e = Q/e- 10~ 2 

\ \ ■— x/e = H/e - 10- 1 

V \ Eq. (38) 
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FIG. 6: Dimensionless specific heat Cg/fcs as a function of 
rescaled temperature ksT/e = l/(/3e) for fi/s = 1/3 and 
various small values of |x/e| — 



V. CONCLUSIONS 

We illustrated the validity of the Jarzynski equality 
and the work fluctuation theorem in the strong coupling 
regime, for the model Hamiltonian (J7J). The central role 
is played by the thermodynamic partition function of the 
open system, that incorporates the interaction of the sys- 
tem of interest with its environment. The influence of the 
interaction is of major importance even in the seemingly 
trivial case in which the system bath interaction Hamilto- 
nian commutes with both the bath and the system Hamil- 
tonians, notwithstanding claims to the contrary [H]. We 
computed the Hamiltonian of mean force for this model 
explicitly and studied its equilibrium thermodynamics. 
In particular we discussed its entropy and its specific 
heat as functions of temperature and other system pa- 



rameters. Like for other strongly coupled system [14L |X5I] 
these two quantities can become negative at low tem- 
perature. Despite entropy and specific heat may become 
negative, they vanish at zero temperature, in accordance 
with the third law of thermodynamics. The only excep- 
tion to this, is for the special value of coupling strength 
X exactly equal to —s, for which the zero temperature 
entropy is equal to fc^ln2. This result is, however, not 
in contradiction with the third law but rather corrobo- 
rates this law; this is so because the two level fiuctuator 
becomes degenerate in this case, as is clearly indicated 
by the Hamiltonian of mean force. 

Interestingly, recent experiments in circuit cavity 
quantum electrodynamics [3l| use a parameter regime 
where negative entropy and specific heat may appear. 
For the architecture presented in [3l[ , these are expected 
below ~ 22 mK and ~ 20 mK respectively. In cavity 
quantum electrodynamics the Hamiltonian in Eq. (1241) is 
obtained from the time-independent Jaynes-Cummings 
model Hamiltonian in the rotating wave approximation 
and dispersive regime [36j . These conditions imply weak 
coupling \x\ <C fi,E, which in fact is the case for Ref. [3l| . 
Whether a Hamiltonian of the type in Eq. ([7]) , with time- 
dependent e(t), and/or possibly strong coupling can be 
implemented with superconducting circuits remains an 
open problem. 



Acknowledgements 

The authors wish to thank David Zueco for fruitful 
discussions. Financial support by the German Excel- 
lence Initiative via the Nanosystems Initiative Munich 
(NIM) and the Volkswagen Foundation (project 1/80424) 
is gratefully acknowledged. 



[1] C. Jarzynski, Phys. Rev. Lett. 78, 2690 (1997). 
[2] G. E. Crooks, Phys Rev. E 60, 2721(1999). 
[3] H. Tasaki, preprint [arXivxond-mat /0009244| (2000) . 
[4] S. Mukamel, Phys. Rev. Lett. 90, 170604 (2003). 
[5] W. De Roeck and C. Maes, Phys. Rev. E 69, 026115 
(2004). 

[6] M. Esposito and S. Mukamel, Phys. Rev. E 73, 046129 
(2006). 

[7] P. Talkner, E. Lutz, and P. Hanggi, Phys. Rev. E 75, 

050102(R) (2007). 
[8] P. Talkner and P. Hanggi, J. Phys. A 40, F569 (2007). 
[9] D. Andrieux and P. Gaspard, Phys. Rev. Lett. 100, 
230404 (2008). 
[10] G. E Crooks, J. Stat. Mech. P10023 (2008). 
[11] P. Talkner, M. Campisi, and P. Hanggi, J. Stat. Mech. 

P02025 (2009). 
[12] C. Jarzynski, J. Stat. Mech. P09005 (2004). 
[13] M. Campisi, P. Talkner, and P. Hanggi, Phys. Rev. Lett. 
102, 210401 (2009). 



[14] P. Hanggi, G.-L. Ingold, and P. Talkner, New J. Phys. 

10, 115008 (2008). 
[15] G.-L. Ingold, P. Hanggi, and P. Talkner, Phys. Rev. E 

79, 061105 (2009). 
[16] G. W. Ford, J. T. Lewis, and R. F. O'Connell, Phys. Rev. 

Lett. 55, 2273 (1985). 
[17] C. Horhammer and H. Biittner, J. Stat. Phys. 133, 

1161 (2008). 

[18] P. Hanggi and G. L. Ingold, Chaos 15, 026105 (2005). 
[19] P. Hanggi and G.-L. Ingold, Acta Phys. Pol. B 37, 1537 
(2006). 

[20] G.-L. Ingold, Path Integrals and Their Application to 
Dissipative Quantum Systems, Lect. Notes Phys. 611, 
1 (2002). 

[21] T. M. Nieuwenhuizen and A. E. Allahverdyan, Phys. 
Rev. E 66, 036102 (2002). 

[22] T. Dittrich, P. Hanggi, G.-L. Ingold, B. Kramer, 
G. Schon, and W. Zwerger, Quantum Transport and Dis- 
sipation (Wiley- VCH, Weinheim, 1998). 



8 



[23] H. Grabert, P. Schramm, and G.-L. Ingold, Phys. Rep. 

168, 115 (1988). 
[24] H. Grabert, U. Weiss, and P. Talkner, Z. Phys. B 55, 87 

(1984). 

[25] A. O. Caldeira and A. J. Leggett, Ann. Phys. (N.Y.) 149, 
374 (1983). 

[26] R. P. Feynman, Statistical Mechanics (Addison Wesley, 

Redwood City, CA, 1972). 
[27] R. P. Feynman and F. L. Vernon, Jr., Ann. Phys. (N.Y.) 

24, 118 (1963). 

[28] J. Teifel and G. Mahler, Phys. Rev. E 76, 051126 (2007). 
[29] P. Hanggi, P. Talkner, and M. Borkovec, Rev. Mod. Phys. 
62, 251 (1990). 

[30] B. Roux and T. Simonson, Biophys. Chem. 78, 1 (1999). 

[31] D. Schuster, et. ai, Nature 445, 515 (2007). 

[32] Max Planck formulated the 3-rd law of thermodynamics 
in rather general terms in his 3-rd edition of (M. Planck, 
Vorlesungen iiber Thermodynamik; Von Veit, Leipzig, 
1910); stating that "the entropy S of any physical body 
reaches a constant value independent of its physical state 
and chemical composition,...., which can be set to S = 
without loss of generality." He naturally at that time 
could not foresee the dependence of entropy at absolute 



zero temperature on symmetry which may imply a de- 
generacy g(N) of the groundstate of a N-particle system. 
This is the case, for example, with spin degeneracy. Such 
symmetry then yields a positive residual entropy at ab- 
solute zero, S(N) = fcsln<;(V), with degeneracy g(N) 
growing no faster than exponentially with N. This in- 
sight came soon after the spin has been introduced into 
the theory of quantum mechanics by Pauli (1924, "Pauli- 
Eq." 1927) and Dirac (1928): Seemingly, this positive 
residual value of entropy at absolute zero temperature in 
terms of symmetry-degeneracy has been introduced first 
by William F. Giauque [33L I34I]. testing the third law for 
hydrogen and ice, note also the credit given to Giauque 
by Linus Pauling [35| . 
[33] W.F. Giauque and H.L. Johnston, J. Am. Chem. Soc. 
50, 3221 (1928); W.F. Giauque, J. Am. Chem. Soc. 52, 
4816 (1930). 

[34] W.F. Giauque and M.F. Ashley, Phys. Rev. 43, 81 
(1933). 

[35] L. Pauling, J. Am. Chem. Soc. 57, 2680 (1935). 
[36] W. P. Schleich, Quantum Optics in Phase Space (Wiley- 
VHC, Berlin, 2001) 



